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PREFACE 


In  this  volume  the  problem  of  visual  detection  of  (low-flying 
aircraft  is  considered  as  part  of  a  general  inve sligatiop.  into  the 
survivability  of  surveillance  aircraft  penetrating  enemy  territory  at 
low  altitude. 

The  work  reported  upon  herein  has  been  sponsored  by  the 
U.  S.  Army  Transportation  Research  Command  (USATRECOM) 
under  Task  Number  9R  38-10-005-01.  It  was  conducted  at  Canadair 
Ltd.  ,  under  Contract  Number  DA  44-177-TC - 793  with  Ool.  W.  F. 
Molesky,  USATRECOM,  acting  as  Project  Officer.  The 
principal  contributor  to  this  phase  of  the  study  was  Mr.  N. 

Kurdyla. 
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SYMBOLS 

Symbol 

A*k  is  the  bottom  projected  area  of  target  in  square  feet 

Af  .  is  the  frontal  projected  area  of  target  in  square  feet 

Ap  is  the  presented  area  of  target  in  square  feet 

As  is  the  side  projected  area  of  target  in  square  feet 

C  is  the  contrast  of  the  target  relative  to  the  background 

at  the  lens  of  the  eye  in  percent 

Cl  is  the  resultant  contrast  which  takes  into  account  the 

effective  contrast  and  the  variation  in  threshold 
contrast  of  an  individual  in  percent 

Ceff=(Al/I)eff  is  the  effective  contrast  when  light  quanta  falls  within 
the  integration  time  and  integration  distance  of  the 
retina  of  the  eye  in  percent 

C0  is  the  intrinsic  contrast  between  target  and  background 

in  percent 

Ct  =  (Al/l)t  is  the  threshold  contrast  of  an  individuals  eye  in 

percent 

d  is  the  instantaneous  minimum  distance  from  the 

target  to  the  line  of  sight  in  miles 

d^j  is  the  minimum  distance  from  the  target  to  the  line  of 

sight  of  the  ifcb  glimpse  and  j**1  flight  path  in  miles 

dij  is  the  minimum  distance  from  the  target  to  the  line  of 

sight  at  the  beginning  of  the  i**1  scan  period  during 
jth  flight  path  in  miles 

djji  is  the  minimum  distance  of  target* s  position  to  the  line 

of  fixation  when  it  intercepts  the  lobe  boundary 
^50  percent  lobe)  in  miles 
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Symbol 


Du 


h 

I 

I' 


A I 


A  I' 

K 

M 

Meff 


n 


Nk 


is  the  minimum  distance  of  target's  position  to  the 
line  of  fixation  when  it  leaves  the  bounds  of  the 
detection  lobe  in  miles 

is  the  total  distance  the  target  traverses  while  within 
the  detection  lobe  in  miles 

is  the  height  of  the  target  above  ground  level  in  feet 

is  the  background  light  intensity  in  Lamberts 

is  the  light  intensity  when  a  target  obscures  part  of 
the  background  in  Lamberts 

is  the  difference  in  light  intensity  reaching  the  eye 
when  only  the  background  light  falls  on  the  eye,  as 
compared  to  the  background  light  falling  onto  the 
eye  with  target  present  in  Lamberts 

is  the  effective  differential  intensity  of  a  flashing 
light  in  Lamberts 

is  the  number  of  sections  around  the  retinal  image 
perimeter 

is  the  total  number  of  quanta  impinging  on  the  retina 

is  the  number  of  light  quanta  impinging  on  the  retina 
that  fall  within  the  integration  capability  of  the 
retina  (when  source  of  quanta  is  moving) 

is  the  number  of  light  quanta  above  the  background 
intensity  required  for  the  perception  of  light  if 
light  quanta  are  absorbed  within  S  and?iof  each  other 

is  the  number  of  points  in  the  kth  section  of  histogram 


TC  -793 


SYMBOLS 


Symbol 

P 


is  the  probability  of  absorption  of  a  quanta  per  trial 


=  P(D/LS)  is  the  probability  of  visual  detection  given  line  of 
sight  to  the  target 


P 


d 


k 


is  the  probability  of  detection  for  the  k**1  interval 
in  the  histogram 


Pn  is  the  probability  that  at  least  n  light  quanta  will  be 

absorbed  in  a  given  region  of  the  boundary 

Fsee  is  the  probability  that  a  nerve  pulse  will  be  triggered 

by  light  quanta  falling  on  the  retina 


Pgee  is  the  probability  that  a  nerve  pulse  will  be  triggered 

by  light  quanta  falling  on  the  retina  for  i**1  glimpse 
or  scan  and  jth  flight  path 


p®eekt 


§tk 

is  the  probability  of  seeing  for  the  L  point  in  the 
k^*1  section  of  the  histogram 


P(*V?;R)  is  the  probability  of  detection  as  a  function  o£  range 
and  angular  position  of  the  aircraft  relative  to  the 
line  of  sight 


R 


is  the  range  in  appropriate  units 


Ri 

*i' 


is  the  range  upon  entering  the  detection  lobe  in  miles 

is  the  range  upon  leaving  the  detection  lobe  in  miles 

is  the  maximum  detection  range  in  miles 

is  the  time  coordinate  at  the  commencement  of  i ^ 
glimpse  during  jth  flight  path  in  seconds 


ij 


is  the  time  coordinate  at  the  commencement  of  i*h 
scan  during  j*'*1  flight  path  in  seconds 
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Symbol 

h+i.  j 


Tf 


T«i 

*  m 


v 


V 


x 


oj 


X,  Y 
Xij’  Yij 


is  the  time  coordinate  at  the  commencement  of  i^1 
glimpse  and  jth  flight  path  in  seconds 

is  the  time  period  allotted  for  the  glimpse  action  of 
the  search  for  all  glimpses  in  seconds 

is  the  time  during  which  target  is  effectively  within 
the  perceptibility  region  for  i^h  glance  and  jth 
flight  path  in  seconds 

is  the  time  period  allotted  for  the  scanning  action 
of  the  search  for  all  scans  in  seconds 

is  the  target  velocity  along  a  straight  line  flight  path 
in  feet  per  second 

are  the  velocities  along  the  X  axis  and  the  Y  axis 
in  feet  per  second 

is  the  meteorological  visibility  in  miles 

is  the  distance  from  the  observer  measured  along 
the  X  axis  at  which  the  flight  path  intercepts  the 
X  axis  in  feet 

are  the  instantaneous  coordinates  of  the  target  in  mile 

are  the  coordinates  of  the  target  when  the  observer 
is  commencing  the  i**1  glimpse  during  the  jt*1  flight 
path  in  miles 


Xi+lj,  Yi+l,j 


are  the  position  coordinates  of  the  target  at  the  end 
of  the  scan  period  or  alternatively  the  beginning  of  the 
ith  glimpse  j**1  flight  path  in  miles 
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SYMBOLS 
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€tff 

Q 

d9/jt 
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s 

tyjt 


are  the  coordinates  of  target  at  the  beginning  of 
scanning  phase  of  the  search  in  miles 

is  the  angular  distance  on  the  retina  within  which  the 
integration  probability  of  the  eye  approaches  one  in 
minutes . 

is  the  expected  number  of  successes,  or  absorptions 
of  quanta  by  the  retina 

is  the  expected  average  number  of  absorptions  of 
light  quanta  from  moving  targets  falling  on  the 
retina 

is  the  expected  average  number  of  absorptions  of 
light  quanta  (from  a  moving  target)  which  fall  within 
the  integrating  capacity  of  the  retina 

is  the  instantaneous  angular  orientation  of  the 
observer  in  degrees 

is  the  rate  of  angular  scan  by  the  observer  in 
minutes  per  second 

is  the  angular  orientation  of  the  observer  when 
target  first  comes  within  visual  range 

is  angular  orientation  of  observer  during  the  i*k 
glimpse  j**1  flight  path  in  degrees 

is  the  off-line  of  sight  axis  visual  angle  in  degrees 

is  the  instantaneous  angular  position  of  the  target 
relative  to  the  observer  in  minutes 

is  the  retinal  angular  velocity  of  the  targets  image 
in  minutes  per  second 
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Symbol 

T 

X 

% 


X 


la  the  time  of  exposure  for  visual  detection  in 
seconds 

is  the  time  constant  of  the  eye  within  which  time 
the  integration  probability  of  the  retina  approaches 
one  in  seconds 

is  the  duration  of  the  on-phase  of  the  flashing 
signal  or  alternatively  the  duration  that  the  aircraft 
remains  within  the  perceptibility  region  in  seconds 

is  an  angle  defining  the  slope  of  the  j**1  flight  path  in 
degrees 

is  a  uniformly  random  variable  defined  between 
0  and  1 
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SUMMARY 


for  use  iJthe  ^lume .preeent®  the  derivation  of  a  .mathematical  model 
or  use  in  the  determination  of  probability  of  visual  detection  of  low- 

vSdeTe  •  ThG  TSUlt  15  the  faCt°r  P(D/Ls)i  Probability  of 

xnves titl  ion  of  vr  Sight’  WMch  Can  be  USed  in  general 

aimed  weapons  8UrV1Vabl1^  Ci  —eillance  aircraft  against  visually 


0<=  f  •  ,  mathematlcal  formulation  of  the  situation  is  based  on 

geometrical  concepts  of  a  two-dimensional  field,  the  known 
characteristics  of  the  human  eye,  and  an  assumed  technique  of 
scanning  by  an  observer. 

To  determine  the  probability  of  visual  detection  in  a  given 
proposed-  a  M°nteiCarl°  aPProach,  utilizing  a  digital  computer,  is 
patlP«  ’  1/  resalts’  taken  over  a  large  number  of  random  flight 
paths,  would  yield  a  probability  of  detection  as  a  function  of  initially 
chosen  parameters.  initially 


1  ™  v/  11U1UC 


of  the  unavailability  to  th.  contracts  of  a  digital  compS’e'r  of' =*'”' 
»d  mal"  ■“Ch  50lUli°”S  °bti““bU  »  »»  economical 


An  extensive  literature  survey  was  carried  out  in  conjunction 
with  this  study.  Although  only  References  1  through  6  were  used  directly 
in  this  volume,  the  others  listed  are  closely  connected  with  the 
visual  detection  subject  and  could  prove  useful  in  further  studies. 
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CONCLUSIONS 


The  problem  of  visual  detection  of  a  rapidly  moving  target, 
although  essentially  complex,  is  amenable  to  analysis  of  the  Monte 
Carlo  type.  The  mathematical  formulation  in  the  present  report 
is  offered  as  an  approach  which  considers  all  of  the  apparently 
significant  parameters  relating  to  this  problem. 


RECOMMENDATIONS 


The  mathematical  model  given  in  this  report  has  not  yet  been 
numerically  evaluated  since  limitations  to  contractor^  computing 
equipment  have  precluded  an  economical  solution. 

To  establish  the  validity  of  the  mathematical  model,  it  is 
recommended  that  a  Monte  Carlo  process  be  carried  out  on  a  high¬ 
speed  computer  to  obtain  the  probability  of  detection  as  a  function  of 
range  for  various  parameters  such  as  target  size,  target  contrast, 
meteorological  visibility,  crossing  distance  and  target  velocity  for 
a  given  search  pattern. 

Some  experimental  results  are  available  in  Reference  1. 

The  model  should  be  checked  using  parameters  for  which  experimental 
results  are  available. 
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INTRODUCTION 


Among  the  problems  related  to  the  study  of  survivability  of 
low-flying  high-speed  surveillance  aircraft  is  the  determination  of 
the  probability  of  visual  detection  given  line  of  sight  F(D/LS)  by 
ground  observers.  The  most  general  approach  to  this  problem  would 
involve  the  study  of  the  interrelation  of  terrain  effects  to  the  dynamics 
of  the  observer-target  system.  This  approach  would  not,  however, 
lead  to  a  general  solution  but  would  depend  upon  the  characteristics 
of  the  terrain  in  which  the  observer-target  system  were  located. 
Before  this  problem  can  be  solved,  it  is  essential  that  the  character¬ 
istics  of  the  obse rver -target  system  in  an  idealized  situation  be 
understood.  The  present  report  is  an  attempt  to  provide  such  an 
understanding  for  an  idealized  situation  where  no  obstructions  exist; 
the  earth  is  assumed  to  be  flat,  and  the  line  of  sight  distance  is 
given  in  plan  only. 

A  mathematical  formulation  based  upon  the  geometry  of  the 
engagement  situation,  the  known  characteristics  of  the  human  eye, 
and  an  assumed  scan  procedure  of  scan-fixate -  scan  is  presented. 

A  solution  by  the  use  of  Monte  Carlo  techniques  for  the 
probability  density  distribution  of  sightings,  within  the  mathematical 
framework  presented,  is  proposed. 
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ASSUMPTIONS 


In  this  study,  the  problem  of  visual  detection  is  only  taken  to 
a  point  where  the  increase  or  decrease  of  light  energy  due  to  a 
target  in  the  field  of  view  causes  a  perceptable  alteration  in  the  rate 
of  triggering  of  nerve  impulses  sent  to  the  brain.  The  problems  of 
noticing  these  changes  once  they  arrive  in  the  brain  are  not 
considered.  Thus^the  vigilance  behavior  of  the  observer  is  not 
taken  into  ac count, and  he  is  always  assumed  to  be  in  the  optimum 
state  of  awareness  and  readiness. 

The  assumptions  within  the  limits  of  the  study  are  given 

below: 

1)  Visual  detection  probabilities  are  required  on  the 
assumption  that  no  obstructions  to  vision  occur  except 
those  directly  related  to  light  extinction  by  the 
atmosphere.  The  earth  is  assumed  to  be  flat  as  far 
as  the  eye  can  see. 

2)  Within  the  range  of  visual  bounds  of  detection  as 
determined  by  atmospheric  conditions  at  a  given  time, 
the  aircraft  is  assumed  to  follow  straight-line 
paths. 

3)  Since  only  low-level  high-speed  missions  are  being 
considered,  the  aircraft  spend  a  large  fraction  of  the 
time  very  close  to  the  horizon  altitude.  Therefore, 

a  two-dimensional  approach  to  the  detection  problem 
seems  justifiable  for  a  good  first-order  approximation. 

4)  To  simulate  a  possible  field  situation,  the  aircraft  is 
assumed  to  enter  the  visual  detection  area  in  a  random 
fashion  both  as  to  orientation  and  as  to  crossing 
distance. 

5)  The  observer  is  assumecj.  to  be  searching  in  a  scan- 
fixate- scan  fashion,  a  fixed  time  being  given  for 
each  operation.  Different  functional  behaviors  can  be 
given  for  the  scan-fixate  portion  of  the  search. 
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An  aircraft  is  assumed  to  be  a  point  source  in  relation 
to  the  geometric  situation.  Where  the  area  plays  an 
important  part  in  empirical  equation,  it  has  been 
included.  No  matter  what  the  orientation  of  the 
aircraft  relative  to  the  observer,  the  linear  dimension 
of  the  target  will  be  obtained  by  assuming  that  the 
presented  area  is  circular. 

While  the  eye  is  fixated  in  given  direction,  the  aircraft 
will  be  represented  by  a  line  element  whose  length 
represents  the  distance  traversed  by  aircraft  during 
the  fixated  glimpse  time. 

During  the  scan  portion  of  the  search,  the  angular 
velocity  of  the  aircraft  across  the  retina  will  be 
assumed  to  be  due  to  the  evefs  scan  rate^ 

Attraction  of  the  observer  by  acoustical  energy  is 
considered  negligible  since  the  observer  is  already 
assumed  to  be  in  state  of  awareness  towards  air¬ 
borne  targets,  and  since  acoustical  energy  is  highly 
nondirectional,  it  serves  only  to  keep  the  observers 
vigilant. 

The  velocity  of  the  aircraft  is  assumed  to  be  a 
constant  for  given  flight  path. 


An  observer,  who  is  visually  searching  for  low-level  hig. 
aircraft,  is  located  at  coordinate  position  (0,0)  as  shown  in  Figur 
All  the  observer's  scan-fixate  angular  positions  and  all  aircraft  a: 
positions  are  measured  relative  to  the  observer  in  a  counterclock 
direction  with  the  positive  X  axis  as  reference  direction. 


i^Glimpse 


Flight  Path 
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Since  many  flight  paths  need  to  be  considered  to  obtain  a 
probability  density  of  detection  as  a  function  of  range,  R,  and  time, 

't  »  let  the  flight  path  shown  in  Figure  1  be  the  When  the  aircraft 

comes  within  a  radius  Rj^  (the  visual  detection  boundary  to  be  defined 
later)  of  observer,  let  the  observer  be  in  the  process  of  glancing  in 
a  random  direction,  0oj.  During  the  progress  of  the  aircraft  along  its 
flight  path,  the  observer  is  simultaneously  going  through  a  set  search 
pattern  of  scan-glimpse  actions.  Thus,  when  the  aircraft  is  located 
at  position  (X^j,  Y^j)  at  time  tjj,  an  aircraft  on  its  flight  path  has 
reached  a  point  in  space  where  the  observer  is  commencing  the  i*h 
glimpse  for  the  still  undetected  target.  The  exact  number  of  glimpses 
will  depend  upon  the  meterological  conditions  which  determine  R m  and 
upon  aircraft  velocity,  v. 

By  the  time  the  i glimpse  is  completed,  the  aircraft  previous¬ 
ly  at  point  (Xy,  Yy)  has  moved  at  constant  velocity  to  the  point 
(X^j ,  Y jj).  At  this  time,  ^ij>  the  scan  operation  comes  into  effect  during 
which  time  the  aircraft  moves  continuously  to  a  new  position.  At  the 
end  of  the  scan  period  the  aircraft  is  located  at  jt  Y^-f^  j)# 

and  the  time  becomes  ^i-fl,  j»  whereupon*  a  glimpse  is  taken  again. 

When  the  observer  takes  the  i**1  glimpse  for  the  j**1  flight  path, 
the  angular  orientation  of  observer  is  denoted  by  ©£ j .  For  i  +  1  s* 
glimpse,  angular  orientation  is  denoted  by  0  •  Similarly,  the 

aircraft1  s  angular  positions  relative  to  the  observer  are  denoted  by 
for  the  beginning  of  the  it^1  fixation  period,  and&^ij  for  the  beginning  of 
the  i*h  scan  period.  With  this  notation  there  is  also  the  freedom  of 
using  general  variables  without  subscripts  to  denote  instantaneous 
positions  of  glimpses,  scans,  and  aircraft  positions  as  is  done  later 
on. 


In  reality,  an  aircraft  flight  path  is  continuous  but  in  this 
discrete  operation  representation  of  the  situation,  the  flight  path  must 
be  segmented  into  line  elements  representing  the  sequential  relations 
between  observer's  search  operations  and  aircraft's  time  positions  as 
illustrated  in  Figure  1,  The  necessity  for  the  separation  of  events  into 
two  discrete  cases  is  necessitated  by  the  characteristic  operation  of 
the  human  eye.  Maximum  acuity  is  obtained  when  the  eye  fixates  a 
given  point  of  space  for  a  period  long  compared  to  a  critical  length  of 
time,  Tc  -  However,  under  the  dynamic  conditions  involved  in  the 
present  study,  the  situation  in  the  sky-target  position  changes 
appreciably  during  the  time  spent  in  re -fixating.  Moreover,  the 
ability  of  the  eye  to  detect  is  not  lost  during  the  scan  portion  of  the 
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search,  but,  instead,  merely  obeys  a  law  strongly  dependent  on  the 
rate  of  angular  movement  of  the  image  on  the  retina.  Thus,  it  is 
deemed  necessary  to  consider  both  effects  for  an  accurate  description 
of  the  search  phases  involved  in  detection.  Also,  by  considering  these 
two  effective  modes  of  detection,  it  may  be  possible  to  specify  an 
optimum  combination  of  glimpse- scan  times  for  greatest  detecting 
efficiency. 

Denote  the  random  flight  path  by  its  angular  inclination  to 
positive  X  axis,  0oj»  anc*  by  its  crossing  point  on  the  X  axis,  xoj. 
During  the  time  of  the  ith  fixation,  Tf,  the  aircraft  moves  a  distance 
vTf  along  its  flight  path  from  the  point  (Xjj,  Yij)  to  the  point 
(X^j-vTf  cos  $oj>  Yjj-vTf  s*n  ♦oj)-  Subsequently,  when  scanning  takes 
place,  the  aircraft  moves  along  from  position  (X^j-vTf  cos  $0j, 

Yij-vTf  sin  ^oj)  to  the  position  [Xij-v(Tf  +  Tm)  cos  (f>0j,  Yij-v(Tf  +  Tm) 
sin  $>oj]  where  Tm  is  the  time  period  allotted  for  the  ith  scanning 
action. 


From  the  above  description,  the  following  relationships 

follow: 


^i+lfj  ~  -^ij  “  v  *  ^m)  cos  ^oj 
Yi+1 ,  j  =  Yij  -  v  (  Tf  +  Tm)  sin  <t>0j 


(1) 


During  the  scan  time,  Tm,  the  observer  sweeps  from  a 
fixation  angle  0£j,  to  the  next  fixation  point  0  j  .  Several 
alternatives  are  available  in  describing  the  behavior  of  a  human 
observer  when  searching.  One  is  to  assume  a  random  glimpse 
between  the  boundaries  of  the  search  area  as  has  been  done  in  most 
previous  reports  on  this  subject.  The  second  way  of  obtaining  a 
description  of  the  re-fixation  process  is  to  assume  that  search 
glimpses  occur  in  a  pattern  transferring  from  one  to  the  successive 
one  by  means  of  a  random  function  which  weights  the  glimpses 
towards  smaller  values.  This  method  or  one  similar  to  it  helps 
simulate  the  realistic  situation  where  an  observer  is  much  more 
likely  to  scan  areas  in  the  vicinity  of  previous  glance  rather  than 
'look1  at  widely  displaced  angles,  unless  other  attention-getting 
information  is  present.  Therefore,  glances  should  be  weighted 
towards  smaller  angular  displacements  and  this  will  be  done  as 
indicated  above. 
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Within  each  segment,  however,  the  angular  position  of  glance  will  be 
determined  by  a  random  variable,  or  a  function  of  a  random  variable, 
to  simulate  the  indeterminacy  of  the  exact  position  of  fixation  points 
when  visually  searching. 

The  main  problem  to  be  solved  here  is  the  probability  of 
visual  detection  of  a  moving  target.  Thus,  it  is  required  to  determine 
whether  or  not  the  moving  point  source  comes  within  the  detection 
visual  angle  (or  field  of  view)  during  the  fixation  period,  and 

«rthermore,  to  ascertain  whether  the  point  target  is  seen  during 
e  time  it  lies  within  the  visual  detection  lobe  limits  while  scanning. 
The  probability  of  seeing,  Pgee#  once  the  target  is  within  the 
observable  off-axis  visual  angle,  depends  on  the  effective  contrast, 
Ceff*  which  is  a  direct  function  of  the  angular  speed  of  the  aircraft 
relative  to  the  observer  (if  a  fixed  time  of  observation  of  target  is 
assumed).  A  similar  situation  exists  while  the  eye  is  scanning  except 
that  here,  the  major  portion  of  the  angular  motion  of  a  target  on  the 
retina  is  a  result  of  the  eyeball  motion.  The  assumption  is  then  made 
that  no  detection  qualities  of  the  eye  are  affected  except  those 
attributable  to  the  rate  effect  of  tho  image  on  the  retina  due  to  eye 
motion  itself. 

To  elaborate  on  the  method  used  to  tackle  this  dynamic 
detection  problem,  reference  must  be  made  again  to  Figure  1.  When 
the  aircraft  is  located  at  point  (Xij,  Yij),  the  distance  from  aircraft 
perpendicular  to  line  of  fixation  of  the  iih  glimpse  is  dij.  At  the 
beginning  of  scan  period,  the  separation  distance  is  d'ij.  At  any 
point  in  between,  let  the  separation  distance  between  aircraft  and 
line  of  ilk  glimpse  be  d.  It  is  now  required  to  compare  whether  the 
distance  subtended  by  the  off-axis  visual  angle.,  £  ,  at  range,  R,  is 
greater  than  separation  distance,  d.  If  so,  there  is  a  chance  for 
detection.  If  not,  no  possibility  for  visual  attention-getting  is 
possible.  Upon  satisfying  the  above  criteria  for  a  chance  of  attention- 
getting,  then  the  probability  of  seeing,  Psee*  which  is  a  function  of 
the  effective  contrast,  Ceff,  and  the  threshold  contrast  Ct  must  be 
included.  However,  the  effective  contrast,  Ceff»  depends  strongly 
on  the  speed  of  the  aircraft,  and  thus  the  detectability  becomes  a 
function  of  the  dynamic  changes  oc curing  during  the  search,  Ceff 
can  also  be  written  as  (  A  I/I  )eff  where  I  is  the  background 
intensity  and  I  =  II  is  the  increase  or  decrease  of  light  intensity 
when  the  target  obscures  part  of  the  background  light. 
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A  THEORETICAL.  MODEL  FOR  DETERMINING  PROBABILITY 
_  OF  VISUAL  DETECTION _ 

■ 

General  Considerations 


A  statistical  probability  detection  distribution  as  a  function  of 
range,  R,  and  of  time  of  exposure,^,  for  each  of  the  parameters  is 
required.  This  can  be  determined  by  using  a  Monte  Carlo  approach. 

To  formulate  the  criterion  of  visual  detection  for  a  Monte 
Carlo  approach,  the  problem  must  be  stated  in  mathematical  terms. 
This  is  done  below: 

The  equation  of  the  jth  flight  path  is 

Y  =  (X  -  xoj)  tan  0oj  (2) 

Similarly  the  equation  of  the  ith  glimpse  becomes 

Y  =  X  tan  9i  j  (3) 

Let  the  distance  from  point  {Xij,  Y i j)  perpendicular  to  the 
line  of  fixation  be  djj. 


Then  , 

dij 


Xij  tan  Qi j  -  Yij 
(1  +  tan2  ©ij)^2 


(4) 


Also,  the  perpendicular  distance  from  point  (Xjj  -  vTfeos  $0j, 
Yij  -  vTfeos  ^oj)  to  the  ith  line  of  fixation,  d 1 1 j  is 

(Xij  -  vTfeos  $oj)tan9ij  -  (Yij  -  vTfsin  $oj) 
d’ij  =  - (1  +  tanZOijjTTZ - - -  (5) 


In  general,  the  instantaneous  distance  of  aircraft  from  line  of 
fixation  is 


Xtan9ij  -  Y 

(1  +  tan20ij)1/2  (6) 

Let  there  now  be  defined  a  quantity  called  the  off-axis  visual 
angle  denoted  by  S  .  This  is  the  maximum  angular  displacement  of 
the  target's  position  from  the  line  of  fixation  within  which  visual 
perception  can  be  considered  greater  than  50  percent.  This  off-axis 
angle  as  show:  in  Figure  2  is  a  function  of  several  parameters  and 
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changes  quite  rapidly  with  range,  R. 


FIGURE  2.  OFF -AXIS  VISUAL  ANGLE 


It  can  be  written  that 

R  tan£>  d  (7) 

is  the  condition  that  must  be  satisfied  for  a  finite  chance  of  detection 
to  occur. 

Also  R£b(X2+Y2)i  (8) 

To  utilize  the  criterion  given  by  equation  (7),  the  functional 
variation  of  the  off-axis  visual  angle,  g  ,  is  required.  A  semi- 
mpirical  formula  has  been  developed  which  combines  the  bulk 

effects  of 

1  light  transmissions  in  the  atmosphere 

2  reception  of  light  by  the  retina 

and  uses  these  to  define  an  off-axis  perceptability  region  as  a 
function  of  the  physical  parameters. 
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To  the  off-axis  visual  angle,  §  ,  are, related  the  basic  factors 
of  presented  area,  range,  intrinsic  contrast  and  meteorological 
visibility  (Reference  4)Ji.  ev 


C0e'3,44R/V=  1.  75  §  *  +  45.  6  R2g/Ap 


where 

C0  is  the  intrinsic  target  background  contrast  in  percent 
R  is  the  range  in  miles 
V  is  meteorological  visibility  in  miles 
§  is  the  off-axis  visual  angle  in  degrees 
Ap  is  the  presented  area  of  target  in  square  feet 

Solving  the  equation  for  £  ,  we  get 


=  < 


3.  0625  +  182.  4  C0(B2/Ap)e 


•3.44B/V12  , 

-i. 


91.  2  R2/At 


"1 


(9) 


(10) 


Equation  (10)  indicates  that  for  a  given  set  of  physical 
parameters  C0>  Ap,  and  V,  the  off-axis  visual  angle,  ,  is  a 
function  of  range,  R.  Thus,  the  value  of  f  at  a  given  R,  defines 
two  points  in  the  visual  detection  lobe.  Repeating  this  for  all 
ranges,  R,  we  then  define  a  visual  detection  lobe  (Figure  3).  All 
points  within  the  detection  lobe  are  assumed  to  have  an  equiprobable 
chance  of  attention- catching  vtlue  equal  to  one  (see  Appendix). 


Thus,  even  though  acuity  changes  as  one  leaves  the  line  of 
fixation,  the  off-axis  visual  angle,  §  ,  defines  the  angular  limits 
where  probability  of  perception  is  reasonably  good.  This  is  so 
everywhere  v/ithin  the  limits  of  the  detection  lobe.  This  effectively 
takes  into  account  the  less  sensitive  parafoveal  and  peripheral 
regions  of  the  retina  as  the  target  comes  within  its  detection 
capability. 

In  order  to  know  where  the  termination  of  the  visual  detection 
zone  occurs,  another  relationship  has  been  presented  by  the  O.  R . 
Evaluation  Group  of  the  Office  of  Naval  Research  (Reference  2). 

Referring  to  the  geometric  situation  (Figure  3),  the  maximum 
range  of  detection,  Rm,  can  be  expressed  again  in  terms  of  the 
physical  meteorological  variables:  i.  e.  , 
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FIGURE  3.  VISUAL  DETECTION  LOBE 


and  where  aLl  others  are  as  previouslv  defined. 


Line  of 
Fixation 


whe  re 


(0.  164)2  Ap  Cc  e'3,  44Rm/v-l.  5?] 
_  Us  ^  ^  Ab  h  +  A#*  (Rm^  - d^  — 

L 


R-m 

Af 

As 

Ab 

h 


is  the  maximum  detection  range  in  miles 
is  the  frontal  projected  area  of  aircraft  in  square  feet 
is  the  side  projected  area  of  aircraft  in  square  feet 
is  the  bottom  projected  area  of  aircraft  in  square  feet 
is  the  height  of  aircraft  above  ground  level  in  feet 
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Thus,  the  distant  boundary  for  visual  detection  varies  with 
the  meteorological  conditions  and  with  the  visual  properties  of  the 
target.  The  direction  of  glimpse  at  the  time  the  aircraft  crosses  the 
boundary  at  Rm  will  be  taken  as  the  reference  point  for  the  definition 
of  a  random  variable,  0oj,  used  in  defining  an  initial  direction  of 
observer's  search  glimpse  for  the  jth.  flight  path. 

Using  the  formulas  presented  up  to  this  point,  the  entire 
geometric  situation  is  defined  for  a  set  position  of  the  i**1  glimpse. 

It  can  be  determined  whether  the  target  satisfies  the  condition  defined 
by  equation  (7)  at  any  time.  If  it  does  satisfy  the  criterion,  then  the 
next  step  is  to  find  the  probability  of  seeing  the  target,  Psee,  given 
that  it  lies  within  the  off-axis  visual  region  given  by  equation  (7). 

Probability  of  Detection 

Once  a  visual  detection  lobe  has  been  defined  by  the  following 

equation 


d  =  R  tan<;  (13) 

it  is  next  required  to  assign  a  probability  of  seeing,  Pgee,  value  for 
a  target  within  the  detection  lobe.  This  step  is  necessary  because 
many  times  in  experimental  situations,  an  observer  will  be  looking 
right  at  the  target  and  will  fail  to  see  it. 

Lamar,  et  al.  in  Reference  3  conducted  a  comprehensive 
study  on  the  detection  of  targets  of  different  sizes,  shapes,  and  contrast 
under  laboratory  conditions.  In  this  paper  they  put  forth  their  theory 
describing  all  their  results  on  a  quantum  statistical  footing.  It,  by  and 
large,  predicts  experimental  results  very  well. 

Included  in  this  theory  is  the  derivation  of  the  probability  of 
seeing,  Psee,  as  a  function  of  the  contrast,  &  I/I.  A  brief  summary  of 
relevant  theory  and  formulas  leading  up  to  this  relationship  is  given 
in  the  following  sections. 

To  explain  the  behavior  of  the  human  eye  On  a  quantum 
statistical  basis,  the  basic  assumption  used  has  been  that  in  order  for 
a  target  to  become  visible,  a  certain  number  of  light  quanta,  n,  must 
be  absorbed  in  some  manner  and  in  some  location  by  the  cones. 
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To  be  more  specific,  it  is  assumed  that  the  area  of  retina  within  which 
n  or  more  quanta  must  be  absorbed  is  further  restricted  to  one  or  more 
sections  of  the  narrow  strip  just  inside  the  retinal  image  boundary* 

What  this  means  is  that  the  presence  of  a  target  will  be  recognized 
if  one  or  more  small  segments  of  the  perimeter  each  absorbs  n  or 
more  light  quanta. 

To  derive  the  probability  of  seeing,  P  ,  from  the  above 
hypothesis,  the  probability  that  at  least  n  light  quanta  out  of  the  M 
light  quanta  impinging  on  the  retina  will  be  absorbed  in  a  given 
section  of  the  boundary  -  i.  e.  ,  Pn  -  must  be  found.  The  probability 
that  less  than  n  quanta  will  be  absorbed  in  the  given  section  is  1-Pn. 

If  there  are  K  sections  around  the  retinal  image  perimeter,  then 
the  probability  that  less  than  n  quanta  will  be  absorbed  in  each  section 
is  (l-Pn)K.  Therefore,  the  probability,  Psee,  that  at  least  n  quanta 
will  be  absorbed  in  at  least  one  of  the  sections  is 

psee  =  i  -  u-pn)K  (14> 

To  find  Pn,  the  probability  of  at  least  n  successes  out  of  M 
trials,  use  is  made  of  Poisson's  Law.  If  p  denotes  the  probability 
of  absorption  of  quanta  per  trial  and  M  the  total  number  of  quanta 
impinging  on  the  retina,  then  the  expected  number  of  successes,  or 
absorptions  of  quanta  by  the  retina,  is 

(:  ~  Mp  (15) 

In  terms  of  the  expected  number  of  successes,  £  ,  the 
probability  of  exactly  n  successes  in  M  trials  is  given  by  Poisson's 
Law  as 

Pn  =  (£ne“*)/n.'  (16) 

The  probability  Pn  of  at  least  n  successes  out  of  M  trials  is 

given  by 

M  £ 

Pn=  ^(£ne  )/n!  (17) 

n 
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Referring  to  equation  (15),  for  finite  values  of  ,  if  p  is 
small,  M  must  be  very  large,  so  that  equation  (17)  becomes 

oo 

pn  =  £(ene'e  )/n.'  (18) 

n 

—  —  Substituting  equation  (18)  into  equation  (14), 

Psee  =  1  -  [l-  e-fe  )/nrj  K  (19) 

n 

When  equation  (19)  is  used  to  predict  the  experimental  curve 
Psee  vs  (Al/I)eff,  it  is  found  that  equation  (19)  describes  the 
experimental  curves  best  when 

n  =  4,  the  number  of  quanta  absorbed  for  perception 
of  target,*  and 

K  — 10,  where  K  takes  on  the  concept  of  the  number 
of  cones  on  the  perimeter  of  retinal  image. 

Thus,  equation  (19)  becomes 

toO  “I  1  0 

1-  2>ne"*)/n!J  (2°) 

The  expected  average  number  of  quanta  absorbed  in  the 
retina,  £  ,  is  proportional  to  the  effective  contrast  between  target 
and  background,  (Al/l)ef£,  and  therefore,  Psee  vs  (Al/l)eff  can  be 
obtained  from  equation  (19). 

The  value  chosen  for  K  depends  upon  the  perimeter  of  the 
retinal  image  but  the  variation  in  Psee  is  quite  small  for  large 
variations  in  retinal  perimeters.  A  value  of  K  =  10  describes  the 
average  probability  of  seeing  curve  to  within  about  5  to  10  percent  for 
values  of  retinal  perimeter  of  interest. 

Thus,  it  is  seen  that  equation  (19)  is  only  a  function  of  £ 

where 

€  =  Kj  (Al/I)eff  (21) 

and  Kj  is  a  constant  of  proportionality. 
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Then, 

Psee  =  1  -  { 1  -  5  [K1  (Al/I)eff] n  e  (AI/I)e££^n.y  (22) 

By  substituting  a  point  on  the  average  probability  of  seeing 
curve,  Kq  ,  may  be  determined.  Equation  (22)  then  gives  a  functional 
description  of  Pgee  vs  (Al/l)e££. 

Because  of  individual  differences  in  threshold  levels  of 
perception,  the  position  of  the  Pgee  vs  (Al/X)e££  curve  shifts  along 
the  abscissa,  the  shape,  however,  remaining  the  same.  To 
normalize  the  curve  for  general  use,  (Al/I)e££  -  (Al/I)t  is  substituted 
for  the  abscissa  instead  of  (Al/l)eff. 

The  contrast  threshold,  (Al/l)^,  can  then  be  obtained  for  an 
individual  and  substituted  into  equation  (22). 

1.  0 
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Tog  (Al/I)eff  -  Log  (Al/I)t 

FIGURE  4.  CHANCE  OF  PERCEPTION  VS  CONTRAST 

The  normalized  curve 

Psee  vs  lo®  (AI/I)ef£  -  log  (Al/I)t  (23) 

shown  in  Figure  4  does  not  change  appreciably  with  daylight  to  dusk 
intensities,  perimeter  of  retinal  image,  and  individual  differences.  It 
is  a  statistically  averaged  curve,  good  for  all  daylight  detection 
conditions. 
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The  value  of  contrast  (Al/l)eff  to  be  used  in  equation  (22)  is 
the  effective  contrast  at  the  retina  of  eye.  Therefore,  if  the 
intrinsic  contrast  between  target  and  background  is  CQ,  then  the  value 
of  (Al/I)eff  used  in  equation  (22)  is  related  to  intrinsic  contrast  as 
follows  (assuming  target  not  moving  in  the  sky): 

C  =  (Al/l)eff  =  CQ  e'3-  912  R/V  (24) 

because  of  the  attenuation  due  to  the  atmosphere. 

Hence,  (AI/I)eff  -  (Al/I)t  =  CQ  e~3'  912  R /V  -  (Al/I)t  *  Cj  (25) 

where 

Ci  is  the  resultant  contrast  to  be  used  in  determining 

probability  of  seeing, 

and  where  all  other  symbols  are  as  previously  defined. 

Equation  (22)  can  be  modified  as  follows: 

Psee  =  1  "  [l  -  5  (KlCl)n  e'KlCl/ni]  (26) 

4 

What  has  been  accomplished  up  to  this  point  is  the  formulation 
of  a  criterion  for  the  visual  detection  of  point  targets.  If  the  target 
is  within  the  detection  lobe,  the  probability  of  seeing,  Pgee,  must  be 
determined,  this  depending  upon  the  contrast,  C,  at  the  retina.  On 
the  other  hand,  if  the  target  is  outside  the  bounds  of  the  visual 
detection  lobe,  the  probability  of  seeing,  Pge^,  is  taken  to  be  uniform¬ 
ly  zero. 


The  next  major  step  required  is  to  incorporate  quantitatively 
the  visual  detection  characteristics  of  a  moving  target.  This  involves 
a  rate  of  movement  in  the  retinal  image  of  the  target  and  hence  the 
integrating  capacity  of  the  eye  must  be  considered.  This  dynamic  aspect 
to  lie  target  detection  problem  will  now  be  examined. 

Detection  of  Moving  Targets 

The  extension  of  the  visual  detection  analysis  to  include  the 
dynamic  effects  of  high-speed  aircraft  was  made  possible  upon  finding, 
and  then  examining,  a  theoretical  treatment  of  contrast  thresholds  of 
moving-point  sources. 
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In  Reference  4  it  was  found  that  the  variable  which  is  effectively 
altered  under  dynamic  conditions  was  M,  the  total  number  of  light 
quanta  absorbed  by  the  retina.  The  quantity  M,  as  will  be  shown 
later,  is  directly  proportional  to  the  contrast,  C,  at  the  eye. 

It  was  shown  in  the  previous  section  of  this  report  that  the 
probability  of  seeing,  Psee,  given  that  target  is  within  the 
perceptibility  region,  was  a  universal  function  of  the  resultant 
contrast,  C^.  Keeping  in  mind  that  the  total  energy  from  a  target 
impinging  on  the  retina,  M,  is  proportional  to  the  contrast,  C, 
the  theory  applying  to  stationary  targets  may  be  extended  to  that 
which  included  the  dynamic  variables  influencing  the  visual 
detection  of  moving  targets.  This  can  be  done  by  inserting  the 
changes  due  to  motion  into  the  resultant  contrast,  C^,  which,  in 
turn,  is  the  abscissa  in  the  universal  curve  P8ee  vs  Cp  in 
Figure  4.  In  this  way  the  probability  of  seeing,  Psee,  value  is 
modified  to  include  the  motion  of  the  target. 

Proceeding  with  the  above-mentioned  approach  of  including 
the  dynamic  effects  of  detection  and  the  quantitative  relationships 
for  the  energy  effectively  available  for  absorption  when  the  retinal 
image  has  a  velocity  on  the  retina  may  be  determined.  A  general 
statement  of  the  quantitative  relationship  between  energy  (number 
of  light  quanta)  required  for  perception  of  light  and  the  angular 
velocity  of  the  target  can  be  stated  as  follows: 

For  a  fixed  flash  time  of  target  the  number  of  light  quanta 
from  the  target  impinging  on  the  retina 

M  Constant  when  ^dfc^/dt')  'Te  ^  <S  (a) 

and  (27) 

MA^(d2/Vdt)  when  ^dT^/dt )  ^  (b) 

where 

'Te  is  the  time  constant  of  the  eye  within  which  time  the 

integration  probability  of  the  retina  approaches  one  (in 
seconds) 

£  is  the  distance  on  the  retina  within  which  the  integration 

probability  of  the  retina  approaches  one  (in  minutes). 
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n  is  the  number  of  light  quanta  above  the  background  intensity 

required  for  the  perception  of  light  if  light  quanta  are 
absorbed  within  TJJ  and  &  of  each  other 

j 

is  the  retinal  angular  velocity  of  target's  image  in 
^  minutes/seconds. 

Equation  (27b)applies  if  the  condition  d  "ft/dt  >>  &  holds,  and 
therefore,  the  relationship  of  M  to  dl^/dt  is  actually  a  description 
of  the  asymptote  to  the  exact  curve  when  it  is  near  dT^/dt  =  &  /?g  . 
Thus,  the  asymptote  relating  M  to  dV*/dt  originates  at  dlfi/dt  =&f 
and  rises  with  a  constant  slope  (n-l)/nif  plotted  in  log-log  co¬ 
ordinates. 


d?^/dt  — 

FIGURE  5.  ADDITIONAL  LIGHT  ENERGY  REQUIRED  WITH 

ANGULAR  VELOCITY 


By  utilizing  the  equations  of  the  asymptotes  rather  than  the 
exact  curve,  the  proportionality  constants  can  be  inserted  into 
equation  (27)  as  follows: 

M  .  Meff(| £)^/  (4)^  »  4 


Meff 


when 
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Meff  is  the  number  of  light  quanta  impinging  on  the  retina 
that  fall  within  the  integration  capability  of  the 
retina  when  source  of  quanta  is  moving, 

and  where  all  others  are  as  previously  defined. 

Using  the  relationship  given  in  equation  (15),  the  expected 
average  number  of  absorptions  of  quanta  (or  successes),  £  ,  can 

be  related  to  M,  the  total  number  of  quanta  impinging  on  retina, 
and  p,  the  probability  of  absorption  of  a  given  quantum. 


=  Mp 

(29) 

=  pMeff 

(30) 

where 

Gz  is  the  expected  average  number  of  absorptions  of  light 

quanta  from  a  moving  target  falling  on  the  retina,  and 

£eff  is  the  effective  expected  average  number  of  absorptions 

of  light  quanta  from  a  moving  target  which  fall  within  the 
integrating  capacity  of  the  eye. 

Using  equations  (29)  and  (30),  we  can  convert  equation  (28) 
to  the  following: 

£  =  eff  (dt^/dt/  (  8  /  T"e)^  ^  when  (d2A/dt)  Ze  £> 

(31) 

£  =  6eff  when  (dZA/ dt)  "Ze  ^  & 

Equation  (31)  has  been  formulated  for  a  fixed  time  of  observation, 
but  in  this  problem  the  time  that  the  aircraft  remains  within  the 
detection  lobe  changes.  In  fact,  under  the  high-speed  conditions  of 
target  motion,  the  time  that  the  aircraft  remains  within  the  detection 
lobe  can  be  represented  as  different  durations  of  flash  intensities. 

An  equation  by  Blondel  and  Rey  mentioned  in  Reference  5  gives 
a  relationship  between  a  flash-type  signal  and  a  steady  signal. 
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Since  La  steady  signal  is  one  whose  duration  is  large  compared 
to and  the  experimental  method  which  will  be  utilized  in  assigning 
a  value  for  n  pre- supposes  an  equivalently  steady  signal,  the  equation 
can  be  incorporated  into  the  proposed  model  to  give  dependence  on 
flash  durations  of  target. 

This  equation  can  be  written  as  follows: 

(A  I')  / (A I)  =  (KTf)/(k  +  Tf)  K  =  1 

k  =  0.  21  (32) 


where 

A  I*  is  the  effective  differential  intensity  of  a  flashing  light 

AI  is  the  effective  differential  intensity  of  a  steady  light 

Tf  is  the  duration  on-phase  of  the  flashing  signal  or  the 

duration  that  the  aircraft  remains  within  the 
perceptibility  region. 

Also,  from  equation  (22)  it  is  observed  that  &  ,  the  expected 
average  number  of  absorptions,  is  proportional  to  retinal  illumination 
A>E/E  and  hence  to  contrast  Al/l  at  the  retina: 

i.e.  ,  6  =  K0(AE/E)  =  K^Al/l)  (33) 

where  KQ  and  Kj  are  constants. 

Keeping  the  above  relationships  in  mind  and  utilizing  the  results 
stated  in  equation  (32), 

( e')/(e)  =  (avIdhai/i)  =  (rf)/(o. 21  +  rf)  (since  11  =  i)  (34) 

Utilizing  equations  (31)  and  (33),  6eff,  the  effective  expected 
average  number  of  absorptions,  may  be  related  to  £r,  the  expected 
average  number  of  absorptions,  when  both  image  velocities  and 
flash  durations  are  taken  into  account. 

n-1 

6eff  _  £eff  _/  S/r  \  n  _ _ 

£  <£'  £  "  \d2/Vdt  /  ‘  0.21  +  Tf 
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so  that 

3  eff 


VlZ^/dt/ 


n-1 

n 


0.  21  +  Tf 


0.  21  4=  TV 


when  7^  2>S] 


when  -gj-  Te  <  *>] 


>'2f  >7i  (35) 


Utilizing  equation  (33),  equation  (35)  can  be  converted  into 
the  following  form: 


-r^=(i?7dr; 

Cpff  _  ^ f 


n-1 


Ceff  _  (&  n 


leil 


0. 21  +  Tf 


o.zTirr  whei’  "3T  ' r*  »  & 


when  Te  <  8 


(36) 


The  contrast,  C,  can  be  further  related  to  the  intrinsic 
contrast  between  target  and  background.  This  is  done  by  introducing 
an  exponential  attenuation  factor  because  of  atmosphere  scattering 
and  refraction  as  light  quanta  travel  towards  the  retina. 

-3.  44R  /V 

C  “  C  Qe  (37) 


Having  related  the  effective  contrast,  Cef£,  involving  retinal 
velocity  and  flash  time  of  target  observation  to  the  intrinsic  con¬ 
trast,  COI  a  physical  characteristic  of  an  aircraft  under  typical  sky 
background,  it  remains  to  obtain  realistic  values  for  n,  dZ/’/dt,  Tk  and 
6  under  daylight  adapted  conditions  of  the  eye.  Each  of  the 
quantities  will  be  considered  separately  and  a  value  will  be  assigned 
where  possible. 


Value  for  n 


As  mentioned  earlier  in  this  volume  in  connection  with  the 
derivation  of  the  probability  of  seeing,  Pgee,  the  number  of  light 
quanta,  n,  that  must  be  absorbed  within  teand  6  under  daylight 
adapted  conditions  was  found  to  be  4  for  the  perception  of  a 
differential  light  intensity  (Reference  5). 
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Therefore,  n  =  4  in  the  relationship  established  in  equation 

(37). 

Value  for  *Ze 


In  a  paper  by  M.  A.  Bouman  (Reference  6),  it  was  found 
that  the  time  constant  of  integration  for  the  eye  can  be  taken  as 
*£&  =  0*05  seconds  with  a  high  degree  of  accuracy.  This  value  for  7i 
is  independent  of  the  wavelength  and  even  the  location  on  the  retina 
upon  which  light  impinges. 

Value  for  <5 


The  value  for  6,  the  effective  length  on  retina  within  which  4 
light  quanta  must  be  absorbed  in  order  to  perceive  an  increment  of 
light,  depends  on  the  wavelength  of  light  and  on  retinal  position 
upon  which  light  falls  (see  Reference  2  for  details). 

However,  under  daylight-to-dusk  light  intensities,  Lamar, 
et  al,  in  Reference  3  indicate  that  the  effective  value  for  S  ifl 
approximately  the  cone  diameter  when  target  detection  is  studied 
for  a  daylight  distribution  of  wavelengths;  i.  ev  &  =0.55  minutes 

of  arc. 

However,  it  is  quite  conceivable  that  6  can  vary  substantially 
from  this  value,  depending  on  the  wavelength  or  wavelengths  reflect¬ 
ed  by  target. 

The  variation  of  6  with  retinal  position  is  assumed  negligible 
because  the  main  detecting  area  of  the  retina  is  the  foveal  region, 
and  there,  the  change  is  negligible. 

To  derive  an  approximation  for  the  time  that  the  aircraft 
remains  in  perception  region, 

let  djjj  be  the  distance  of  aircraft  position  perpendicular 
to  line  of  fixation  when  it  intercepts  the  lobe 
boundary 

be  the  distance  of  aircraft  position  perpendicular 
to  line  of  fixation  when  it  leaves  the  bounds  of  the 
detection  lobe 
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D 


be  the  total  distance  aircraft  traverses  while 
within  detection  lobe 


be  the  range  upon  entering  the  detection  lobe 
be  the  range  upon  leaving  the  detection  lobe, 
and  let  all  others  be  as  previously  defined. 


*lf 


D 


ij 


=  (di j i  +  d'ij,>  cos  C^72  "  {0ij  ~  M] 


Also, 


lJl 

di  j  =  R  ^  tan  § 

d'iji=Ri'tan^ 


Dy  -  ^Rjtan§  +  R^1  tan  (•'  j  sin  (0jj  - 


(ft 


£* 


.  .  e;  +  i  «->  a  q  ro  .  ^  /j^  \  e  / 

1^..  .  * - *  '-O'  1  "‘'■p'  c 


91. 2  (R  x  -  /  Ap) 


(r  2  -3.  44R  i  1  /VT 

I  3.  0625  +  182.  4  CQ  (Rj:  /Ap)  e 


91. 2  (Ri'2/Ap) 


-  1,  75 


(38) 

(39) 

(40) 

(41) 

(42) 


Therefore,  the  time  during  which  the  target  is  effectively  within 
the  perceptibility  region  for  i ^  glance  and  j***1  flight  path  is 


Tfij  =  Dy/v  =  (l/v)  (Rj  tan  %1+Ri'tan  )  sin  (0ij  -  <|)0j) 


(43) 


Value  for  d#7dt 


To  apply  equation  (36),  it  is  necessary  to  obtain  dlf/dt  as  a 
function  of  aircraft  velocity  and  aircraft  position  co-ordinates. 

For  the  aircraft’s  instantaneous  position  (X,Y)  as  shown  in 
Figure  7, 
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FIGURE  7.  ANGULAR  VELOCITY  RELATIVE  TO  OBSERVER 


7/ 

*  57.  3  x  60  tan-1  Y/X 

(44) 

dV 

dt 

=  57.  3  x  60 

X *  +  Y^ 

(45) 

vx  and  Vy  are  velocities  along  the  X  and  Y  axis 

d^/dt  is  the  rate  of  movement  of  retinal  image  in  minutes/ 
second . 


Armed  with  discrete  values  of  Z“e,  &  ,  n#  and  function  value  of 
Tf  and  d2fi/dt  together  with  the  foregoing  comments  as  to  the  variations 
of  the  above  parameters,  equation  (36)  becomes: 

,  11  2N3/4 

Ceff  -  C  Qe  ~3- 44^l/yf  57.  3  x  6°  fjl  \  gytan£l  +  Ri'tangl'jsinte  -  0) 

CY  l(R ^tan^i  +  Ri,tan%i')sin(G’  -  0)  +  0.  21v 

7  (46a) 


,vyx  " 


when  dZVdt  >  S/'te  =  H  minute s / second. 
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(R  1  tan  +  R  l  1  tan  5l  *)  sin  (e  -  (ft) _ 

(R  i  tan  +  R  i  1  tan  cos  {©-  (J))  +  O.Zlv  {46  b) 

=  11  minute  s /second. 

Separation  of  Visual  Detection  Area  into  Regions  dZ/Vdt  >  S/2e  jd^/dt<6/Z^ 

To  find  the  geometric  regions  within  which  *2!^  d zH  dt  ^  , 

equation  (45)  is  manipulated  by  substituting  vsin(t>0j  for  vy  and  vcos(l)0j 
for  vx: 

dt^/dt  =  57.  3  x  60  v  (X sin  0oj  -  Ycos  C{>oj)/(X2  +  Y2)  (47) 

Then,  57.  3  x  60  v  T<>  (Xsin  0oj  -  Ycos  <J>oj)/(X^  +  Y2)  ^  £> 

is  the  region  which  must  be  satisfied,  where  S  is  in  minutes. 

57.  3  x  60  v(  f  )  (Xsin  0oj  -  Ycos  Cj>0j)  %  X2  +  Y2 

X2  -  3438  v ( ) X s in  <J>0j  +  [1  "7 1  9  v(^£  )  sin  <t>0j~|  2  +  Y2  + 

3438  v(-|)Ycos  0>oj'  +  J1719  v(-j&)  cos  <{)ojj  2  4  (1719  v-2l)2 

|X  -  1719  v(^)sin  (t)0jj2  +  [Y  +  1719  v(|-)cos  0OJ|  2  <  (1719  v-gt)2 

(48) 

The  locus  described  by  equation  (48)  is  a  circle  with  center  at 
IT" 719  v(^-)sin  Ooj.  -  1719  v(-^)cos  (pjjj  anc*  radius  of  length  1719* 

vz?/a.  -1 


-3.  44R,  /V 

Ceff  =  CDe 

when  dZ^/dt  ^  S I 'te. 


This  locus  gives  the  boundaries  separating  the  area  into  regions 
where  d^/dt  )>  SJ2c  and  d^/dt  <  6>/  TTc.  ,  different  effective  contrasts, 
Ceff,  applying  in  each  of  these  regions  Qsee  equation  (46^. 

Assuming  a  typical  value  of  the  linear  velocity  of  the  aircraft 
(i.  e.  ,  v  =  1000  feet/second  or  1/5  mile  /  second) ,  an  order  of  magnitude 
for  the  radius  of  the  circle  can  be  found. 

Radius  of  circle  =  1719  v  7e/&  =  34.  3  miles 
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Then, 


Center  of  circle  =  34.  3  sin  <t>oj>  -  34.  3  cos  <J>0j  miles 

These  regions  are  drawn  graphically  for  the  jth-  flight  path 
and  i**1  glimpse  in  Figure  8. 


jth  Flight  Path 
Te(dZ//dt)  <  8 


7r_(d^/dt)<& 


Visual 
Bounda  ry 


FIGURE  8.  REGIONS  WHERE  ANGULAR  VELOCITIES 

ARE  GREATER  THAN  AND  LESS  THAN  6/T© 


It  is  seen  that  most  of  the  area  in  which  visual  detection  is 
possible  involves  the  c ondition  (dz^/dt)^  >  &  ,  and  therefore,  equation 
(46a)  can  be  used  uniformly  to  obtain  Ceff. 

All  the  concepts  and  formulas  considered  have  been  directed 
towards  describing  the  visual  detection  ability  of  the  fixated  eye 
looking  for  a  moving  target.  The  same  concepts  can  also  be  applied  to 
the  scanning  search  mode  of  the  eyes.  It  is  assumed  that  the  retinal 
angular  velocity  of  the  target  is  due  entirely  to  the  relative^  angular 
velocity  of  the  observer's  eyes  and  that  no  other  acuity  changes  occur 
in  the  retinal  detection  ability.  *' 

Then  in  equation  (36)  dlf'/dt  can  be  replaced  by  d©/dt  in  degrees, 
the  angular  rate  of  observer  's  scan;  i.  e.  , 

Ceff  =  C-0e  -3*44  R /V  (11/3438  de/dt)  (rp/ 0.21  +Tp)  when 

(d^/dt)r^>  6.  ,0  (49) 
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Here  d©/dt  can  be  chosen  according  to  any  scan  angular 
velocity  function;  j.  e.  ,  a©/dt  »  constant  is  the  time  of  observation  of 
target  (time  within  detection  lobe  limits  during  scanning). 

2  g  =  2  f[3.  0625  +  182,4  C0(R2/Ap)]^  -  1.  75~)  2 

d©7dt  d©7dt  ^  91.  2H2/Ap  j  (50) 

Description  of  Glimpse  and  Scan  Behavior  by  Observer 

From  the  following  two  modes  of  possible  detection, 

1)  the  fixated  glimpse  mode  in  which  the  maximum  acuity 
of  the  eye  is  utilized  (the  achievement  of  maximum 
acuity  is  coupled  with  a  relatively  small  solid  angle  of 
coverage  per  unit  time),  and 

2)  the  scanning  mode  in  which  acuity  is  reduced  (but  on  the 
other  hand,  the  solid  angle  swept  can  be  extremely  large) 
and  which  is  required  to  define  a  scanning -fixate  pattern  for 
target  search, 

a  pattern  of  behavior  of  a  human  observer  under  alerted 

type  conditions  will  be  postulated.  The  main  properties  that  will 

be  included  in  the  scanning-glimpse  pattern  will  be  as  follows: 

1)  The  necessity  of  including  a  randomness  in  the  selection 
of  the  line  of  visual  fixation  of  the  observer. 

2)  The  characteristic  of  most  human  observers  to  glimpse 
at  an  angle  displaced  from  original  glimpse  by  a  most 
probable  angular  displacement. 

A  function  which  is  realistic  and  yet  fairly  simple  can  be 
defined  to  help  position  successive  fixated  lines  of  vision. 

This  can  be  written  on  an  empirical  basis  as  follows: 

©i+l,j  "  ©ij  =  ( 7772 )  sin(2  X  i j  -  1)(1T/Z)  when  0  $©i+1  f  j  <  77" (51) 

where 

^  j  is  the  ar  gular  position  of  glimpse  relative  to  positive 
X  axis  for  i  b  1st  glimpse,  jth  flight  path 
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©ij  is  angular  position  of  glimpse  for  i**1  glimpse  flight  path' 
Xy  is  a  uniformly  random  variable  defined  between  0  and  1. 


FIGURE  9.  POSITION  OF  OBSERVER'S  MOST  LIKELY  GLANCE 


By  choosing  a  random  number  between  0  and  1  forXjj*  we  get 
a  most  likely  position  for  the  successive  glimpse  according  to  a 
sinusoidal  function  {Figure  9). 

To  keep  this  dynamic  analysis  simple,  it  is  further  assumed 
that  the  rate  of  angular  rotation  when  scanning  towards  a  new  fixation 
is  constant;  i.  e.  , 

d©/dt  =  c  (52) 

Hence,  the  time  that  the  taiget  remains  within  the  perception 
region  is  constant  as  the  target  image  moves  across  the  retina, 

3 /4 

i.e.  ,  Cef£  =  (U/c)  •  r£/(0.  21  +r£)  (53) 
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A  quick  calculation  indicates  that  any  value  of  d©/dt 
between  11  minute  s  /  second  up  to  a  maximum  of  1  5  degree  s  /  second  is 
within  detection  possibility  when  scanning  if  the  intrinsic  contrast 
between  target  and  background  is  of  the  order  of  30  percent. 

The  value  of  T  #  time  of  scanning  period,  thus  changes  for 
each  glimpse  change  depending  upon  angular  distance  traversed  by  the 
observer  if  equation  (41)  is  obeyed. 

Any  alternative  assumption  of  observer’s  behavior  can  be 
included  if  more  definite  information  is  obtained. 

Method  To  Be  Used  in  Constructing  vs  R  and  vs  'V  Curves 


By  using  a  Monte  Carlo  approach  on  the  computer  for  random 
flight  paths,  the  number  of  detections,  N,  between  each  range,  R,  and 
R  +  AR  can  be  acc\.imulated  and  stored.  From  this  information,  a 
histogram  of  the  number  of  detections,  N,  may  be  made  between  R  and  R 
4*  AR  for  each  range,  R. 


R 

FIGURE  10.  HISTOGRAM  OF  FREQUENCY  OF  DETECTION 

AS  A  FUNCTION  OF  RANGE 

To  obtain  a  probability  of  visual  detection  curve  as  a  function 
of  range,  R,  each  N  must  be  weighted  with  a  probability  of  seeing 
quantity,  ^see»  because  each  possible  detection  has  a  definite 
probability  of  detection.  For  the  k^h  interval  of  histogram, 
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pdk 

where 

pdk 


is  the  probability  of  detection  for  the  k*k 
interval  in  the  histogram 


(54) 


P8ee  is  the  probability  of  seeing  for  the^*1  point 

in  the  k^1  section  of r  the  histogram 

is  the  number  of  points  in  the  k**1  section  of  the 
histogram 

The  exact  subdivision  required  to  obtain  a  well  defined 
probability  of  detection  curve  P4  vs  R  is  not  known.  However, 
sufficient  resolution  should  be  obtained  if  internals  in  R  are  about 
1, 000  yards.  . 


To  obtain  probability  of  detection^Pd^s  the  time  of  exposure, 
'V  t 't'  must  be  related  to  R. 


Defining  *2T  *  0 


when  R  *  Rm  (see  Figure  11)> 


FIGURE  11.  CONSTRUCTION  TO  DETERMINE  TIME  OF 
EXPOSURE  IN  SEARCH  AREA 
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from  the  triangle  AOB,  (vf)2  =  R2  +  ~  2RRmcos  (©  -  0').  (55) 

If  detection  occurs  at  range,  R,  then  time  of  exposure,  'U , 
can  be  derived  from  equation  (43), 

T  =(l/v)  [r2  +  .  2RKmcos(0  -  0')  ]  7  (56) 


whe  re 

X  is  the  time  of  exposure  of  the  target  from  the  point 
it  enters  the  visual  detection  field  to  point  of 
detection  which  is  equivalent  to  the  time  elapsed  since 
aircraft  entered  detection  field  to  the  time  of  detection 

0  is  the  aircrafts  angular  position  relative  to  the 
observer  in  degrees. 

Substituting  for  0  and  0* 

T2  =(l/v2)  [r2  +  Rm2  -  2RRmcos(tan_1  Y/X  -  tan"1  Ym/Xm)]  (57) 

T  =(l/v  )  [r2  +  Rm2  -  2RRmcos(tan_1  Y/X  -  tan-1  Ym/Xm)j  1  (58) 

where 


^m  x m  * * 

Thus,  every  time  a  finite  probability  of  seeing,  Psee,  occurs 
at  range,  R,  then  by  using  equation  (19)  a  Pgee  ^  corre spending  to  each 
T  may  be  determined. 


The  probability  of  detection,  Pdk*  can  ^e  determined  for  the 
interval  in  histogram  between 'V and  'TT  +  AT  by  the  following: 

Nk 


dk 


I 


iek£ 


/Nk 


(59) 


JL  =  1 

From  the  outlines  of  the  above  model,  the  probability  of 
detection,  P^,  vs  range,  R,  and  also  the  probability  of  detection,  P^, 
vs  time  of  exposure,  ,  as  a  function  of  all  the  parameters  involved 
may  be  obtained. 
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APPLICATION  OF  THEORY  INTO  THE  PROBLEM  OF 
SURVIVABILITY  OF  MANNED  AIRCRAFT 

One  of  the  probability  factors  in  the  overall  probability  of 
kill  is  the  probability  of  detection,  given  line  of  sight,  P(D/LS) 
(see  Volume  IV).  In  this  case  line  of  sight  to  the  target  has  been 
assumed  right  from  the  start. 

The  results  that  can  be  obtained  from  the  histogram,  i.  e.  , 
Pd  vs  R,  are  thus  the  ones  representing  P(D/Lg)  for  the  visual 
detection  of  aircraft.  They  can  be  utilized  in  determining  the 
probability  of  kill,  P(K),  of  an  aircraft  when  using  vision  for 
initial  detection  of  the  target. 
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APPENDIX 

The  formula  given  by  equation  (10)  gives  the  off-axis  visual 
angle,  ^  ,  within  which  50  percent  of  the  detections  occur.  It  is 
thus  an  average  angle  giving  the  field  of  view  for  given  aircraft  and 
meteorological  parameters.  It  has  been  assumed  in  the  main 
body  of  this  report  that  all  points  within  the  off-axis  angle  have 
a  weighting  value  of  unity  and  ail  points  outside  this  angular  region 
a  weighting  value  of  zero. 

A  better  approximation  to  the  exact  value  can  be  obtained  by 
assuming  an  off-axis  weighting  distribution  which  is  Gaussian  with 
the  average  RMS  angle  defined  in  terms  of  ^ 

where 

i 

_  -  3  44R  /V  2 

I  [3.  0625  +  182.  4  C0(R"/Ap)  e  _ )  -  1.  75 

’  L  91.2(R2/Ap) 

which  is  a  repetition  of  equation  (10). 


] 


Then,  the  probability  of  detection  as  a  function  of  the  angular 
position  and  range  of  the  aircraft  relative  to  the  line  of  sight  can  be 
expressed  by  the  equation 


P(^,R)  =  e  -i(Btan2^/Rtan  1.43  5)2 


(60) 


where 

P(t/L,R)  is  the  probability  of  receiving  enough  light  energy 

at  the  retina  of  the  eye  for  perception  as  a  function 
of  off-axis  angle  and  range. 

If  distant  targets  are  considered,  then  %  and7^(units  in  degrees) 
in  general  are  small  and  the  equation  reduces  to 

The  off-axis  inefficiencies,  P(Z^,R),  can  be  determined 
experimentally  to  see  how  accurately  a  Gaussian  distribution  describes 
the  experimental  results.  , 
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This  can  be  done  by  using  a  stationary  target  with  a  fixed 
effective  contrast  to  maintain  the  chance  effect  of  seeing,  Psee, 
constant.  Then  by  many  trials,  at  different  off-axis  positions, 
the  number  of  detections  made  to  the  total  number  of  trials  at  that 
off-axis  position  would  give  the  probability  of  detection  for  that 
angular  position.  This  can  be  repeated  for  all  off-axis  positions 
to  check  how  closely  the  experimental  results  conform  to  the 
equation. 

By  including  the  sophistication  for  off-axis  angle  given  in 
this  appendix,  the  probability  formulation  for  visual  detection 
becomes 

P(D/LS)  =  P  Q>",  dl^/dt, r]  =  P('tf',R)  .  Psee(dV7dt,R)  (62) 

Since  range,  R,  is  very  insensitive  over  one  interval  of  the 
histogram,  if  the  intervals  chosen  are  small  enough,  the  above 
probabilities  in  that  interval  are  approximately  independent  and 
hence  can  be  multiplied  together. 

One  then  obtains  the  wanted  formulation  by  using  equations 
(22),  (10)  and  (61) 


P(D/LS) 


where 


and 


Ceff 


P(^,dl>7dt,  R) 

r. 


,  _  .n  -ki  Ceff 
1  _  (kj  Ceff)  e _ 


-M 


435)' 


Ceff  *  Ceff(d1>/dt,R) 

=  Coe-3'44R/V/ 


3/4 


(63) 


Idlk/dt 


0  .  21  +  Tf 

when  dlT/dt  >  11  minute s  /  second 


C0e  ^«44R/\  # _ £ _  when  d~&! dt  <  11  minutes/ 

0.  21  f  Tf  second 
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